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Abstract 

We present a detailed non-perturbative analysis of the time-evolution of a well- 
known quantum-mechanical system — a particle between potential walls — describing 
the decay of unstable states. For sufficiently high barriers, corresponding to unstable 
particles with large lifetimes, we find an exponential decay for intermediate times, 
turning into an asymptotic power decay. We explicitly compute such power terms in 
time as a function of the coupling in the model. The same behavior is obtained with 
a repulsive as well as with an attractive potential, the latter case not being related 
to any tunnelling effect. 
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1 Introduction and summary of results 



The decay of a particle originates from a coupling of its states, by virtue of some interaction, 
to a continuum of many-particle states, which constitute the decay products. It is necessary 
to have a coupling to an infinite number of states because, with a finite number of them, 
one only obtains an oscillatory behavior of particle amplitude with time, as it occurs for 
example with spin systems [I]. Due to its intrinsic complexity, this process is usually 
treated in quantum field theory, as well as in quantum mechanics, in perturbation theory. 
In lowest order, i.e. in the limit of vanishing coupling, one has a stable particle belonging 
to the spectrum of the free theory. By switching on the interaction to the multi-particle 
continuum, the particle becomes unstable and disappears from the spectrum. A common 
property of unstable particles is the exponential decay law, which can be considered a law 
of nature: 

N(t) = N e-^, (1) 

where N(t) is the number of particles at time t, Nq = N(0) is the initial number of 
particles, and r is the mean lifetime. Violation mechanisms of the exponential decay law 
have been proposed by various authors [3], but have escaped experimental detection up to 
now. Such effects manifest themselves in the form of power terms in the time, presumably 
with a small coefficient. In this note we explicitly calculate such power effects in a simple 
quantum-mechanical model, as a function of the interaction strength. We consider a model 
similar to the one originally proposed to explain the a-decay of nuclei: a particle initially 
confined between large and thick potential walls [I]. Because of tunnelling effect, the 
wavefunction will filter through the potential walls to an infinitely large region, so that the 
probability for the particle to remain in the original region will decrease with time down to 
zero. Even in our simple model, the time evolution cannot be computed in closed analytic 
form and we derive rigorous asymptotic expansions for the wavefunction of the unstable 
state at large times; let us stress that the expansion parameter is \jt and not the coupling 
of the interaction. 

The main findings of our work are the followings. For \g\ <C 1, where g is the coupling, 
we find an exponential decay with a large lifetime r « 1/g 2 in the time range 

l«*<!^, (2) 
T 

which turns into an asymptotic power decay ~ g A /t 3 (for l^l 2 )- By increasing \g\, the ex- 
ponential time region reduces and eventually shrinks to zero. That suggests that violations 
to the exponential decay law should be easier to detect in strong coupling phenomena. It is 
remarkable that a similar decay pattern is found for a repulsive as well as for an attractive 
potential. That was not expected on physical grounds, because there is no tunnelling effect 
for "negative walls" . In other words, tunnelling is not the decay mechanism in our model. 
In general, the decay properties seem to originate from the resonance characteristics of the 
system, which resembles a resonance cavity for small \g\ and occur both in the repulsive 
and in the attractive case. 
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2 The model 



The Hamiltonian operator of our model reads: 

where m is the particle mass, A is a coupling constant and x = L > is the support of 
the potential. We assume that wavefunctions ip(x) are defined on the positive axis only, 
x > 0, and that vanish at the origin: ip(0) = 0. Formulae can be simplified by going 
to a proper dimensionless coordinate via x = Lx'/n and rescaling the Hamiltonian as: 
H = Ti 2 u 2 /(2mL 2 )H' . The new Hamiltonian reads 



and contains the single real parameter g = % 2 /{2mLX). Let us omit primes from now on 
for simplicity's sake. 



3 The spectrum 

The eigenfunctions of the Hamiltonian operator H read: 

Mx) « ( - % - e lkx + % - e" fc ) 9(tt - x) + (a k e lkx + b k e~ lkx ) 6{x - tt) , (5) 

and have eigenvalues e k = k 2 . 6(x) — 1 for x > and otherwise is the step function and 
the coefficients a k and b k have the following expressions: 

In general, A; is a complex quantum number. For real g, (a k )* = b k *, implying that the zeros 
of the equation a k = are complex conjugates of the ones of b k = 0. Since b_ k = —a k , 
the zeroes of a k = are opposite to those of b k = 0. Furthermore, the eigenfunctions 
are odd functions of k, i.e. ip- k (x) = —ip k (x), implying that it is necessary to consider 
only "half" of the complex /c-plane. However, in order to have real eigenvalues, as it 
should for the hermitian operator H, k must be either real or purely imaginary; the former 
case corresponds to non-normalizable eigenstates with positive energies in the continuous 
spectrum, while the latter case corresponds to normalizable eigenfunctions with negative 
energies in the discrete spectrum, if any. Since we have to deal both with normalizable and 
non-normalizable eigenfunctions, normalization will be considered case by case in the next 
sections. 
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3.1 Continuos spectrum 



The continuous spectrum is obtained for real k. In this case b k = (a k )* and all the eigen- 
functions are real. Let us normalize them as: 



f 

Jo 



ip* k ,(x)ip k (x)dx = 5(k-k'), (8) 

where S(q) is the Dirac 5-function. The normalization factor reads: 

1 _ /2~ 1 

k{9) ~ y/27ra k b k ~ V tt ^ 1 + \jiixgk) sin2/br + 1/(2ttV^ 2 )(1 - cos 2/bff ' 

The final expression for the eigenf unctions therefore can be written as: 
1 



(--^-e lkx + -^-e~ lkx ) 9(ir-x) + ( ^ e lkx + ^ e~ ikx ) 9(x - tt) 
V 2|a fc | 2|a fc | / \\a k \ \a k \ J 

(10)" 

Note that, because of continuum normalization, the amplitude of the eigenfunctions outside 
the wall is always 0(1), no matter which values are chosen for k and g, while inside the 
cavity the amplitude has a non-trivial dependence on k and g. Note also that we can 
assume k > 0, implying that there is no energy degeneracy. 



3.2 Discrete spectrum 

Let us now consider the eigenfunctions with a purely imaginary quantum number k, k = ik2 
with &2 real, i.e. with the negative energy e ik2 = —k 2 2 < . It holds: 

AM oc ( - \ e~ k * x + \ e +k * x ) 6(n - x) + (a lk2 + b lk2 e +k * x ) 6{x - tt) . (11) 

In order to obtain a normalizable state, the exponentially growing term for x — > +oo must 
vanish, i.e. have zero coefficient — quantization condition. We can impose for example: 

b ik2 = and k 2 > 0. (12) 

The equation is easy rewritten as: 

e -27rk 2 = l + 2ll k 2 g. (13) 

The equation e~ s = 1 + gs, with s = 2nk 2 , has no solution for g > 0, in agreement with 
physical intuition: there are no bound states with a repulsive potential. There is instead 
one non-trivial solution k 2 (g) for — 1 < g < 0, again in agreement with physical intuition: 

Ak 2 (x) = C k2 [9(n - x) (e k * x - e^ x ) + 9{x - tt) (e 2 ^ 2 - l) e^ 2 *] (14) 

By normalizing to one, the overall coefficient reads: 



For — 1 < g < 0, the trascendental equation for k 2 has the approximate solution 

k 2 (g) ^ 9 -±l + 0[(g + l) 2 ] (16) 

while for a negative coupling of small size, — 1 <C g < 0, 

k 2 (g) = - J-[l-e^ + 0(e 2 ^)l. (17) 

Let us note that by imposing = and k 2 < 0, one obtains the complex-conjugate zero 
—k 2 (g), in agreement with the symmetry property of a k and b k ; the latter zero gives rise 
to the same eigenfunction. 

3.3 Normalizable eigenfunctions with complex energy 

In order to study the temporal evolution of general wavef unctions, it is convenient to 
consider normalizable eigenfunctions with a truly complex k = k± + ik 2 , i.e. with the 
complex energy e = k\ — k 2 +i2kik 2 . The temporal evolution is controlled by the exponential 
factor 

e -iet _ e -i(kf-k%)t+2k 1 k 2 t 

implying decay with time for k\k 2 < 0. The eigenfunctions are of the form: 

i> M2 {x) oc [ - % - e^- k ^ x + % - e (-*+^] e(n - x) + 

+ [a kl+lk2 e^- k ^ x + b kl+tk2 e (-*+ fc 2 )-] 6(x - tt) . (19) 

As in the previous section, normalizable eigenfunctions are obtained by killing the expo- 
nentially growing terms for x — > +oo. By avoiding also an exponential growth with time, 
one obtains the relations: 



'k\+ik2 



= 0, ki < 0, k 2 > 0. (20) 



The above equation has a countable set of solutions k^ n \g) for any real g. For small g (n 
is a positive integer): 

k {n) (g) = -n + gn- g 2 n + ing 2 n 2 + O (g 3 ) . (21) 

The distance of k^ n \g) from the real axis oc g 2 . For small g, the eigenfunctions are confined 
inside the potential barrier, with an exponentially-decaying tail for x > n. 



4 Limiting cases 

The decay properties of our model become simpler in the limiting cases corresponding to 
a large potential barrier, \g\ 1, and to a small potential barrier, \g\ ^> 1. 
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4.1 Low Potential Barrier 

In the free limit, g — > ±00, we have that — > —i/2, bk — > i/2 and the eigenfunctions 
reduce to sinusoidal waves in the whole positive axis: 

[2 

ipk{x) = y-sm(kx) (g = ±00) . (22) 

There is clearly no discrete spectrum in this limit. 

4.2 High Potential Barrier 

Let us consider an eigenfunction in the continuous spectrum in a high potential barrier; 
formally we take 

k « 1 and \g\ < 1 , (23) 

so that 

\g\k < 1. (24) 

For most values of k, 

1 — r S3 |0|A; < 1, (25) 

\ak\ 

i.e. the wavefunction amplitude is much smaller inside the potential wall than outside, 
where it is always 0(1). Let us now determine the values of the quantum number k, if 
any, for which the eigenfunctions have a large amplitude, i.e. an amplitude 3> 1, inside the 
potential barrier (i.e. for < x < it). By writing 

k = n + e (n) , (26) 

with n a non-zero integer, and imposing l/|ajfe| to be as large as possible 0, we obtain: 

k {n \g) = n- gn + O (g 2 ) . (27) 

For these specific values, the amplitude inside the barrier (0 < x < ir) is rather large 
because 

1 „ / 1 



O (-—) > 1 . (28) 

lOfcwi \\g\ n J 

In general, the amplitude is large inside the barrier for values of k inside intervals centered 
around k^ (g) of size 

Sk^{g) ^ \g\n. (29) 

For these ranges of k, the region between the potential walls resembles a resonant cavity. 
On the contrary, for values of k outside the intervals (n—\g\nn, n+\g\n), the wave amplitude 



1 There is no real k which makes bk (or a^) exactly vanishing for g ^ 0. As shown in the previous 
section, the expansion of the solution of the trascendental equation bk = for small g contains an imaginary 
part in second order in g. 

2 For = for example we obtain a pure sinusoidal wave in the whole positive axis, i.e. equal 
amplitudes inside and outside the barrier. For these wavelenghts, there is actually no effect of the potential 
(transparency) . 
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is small inside the barrier, i.e. it is <C 1. Finally, there is a large phase shift when the 
inside amplitude becomes <C 1 from ^> 1 and vice- versa j5]. 

Let us stress that a similar resonant behavior is found for g > as well as for g < 0, as 
it stems from the above formulae. Since 



2tt 



m{g) 



2tt 
n 



l + g + 0(g 2 ) 



(30) 



the only difference is that for g > the "resonant" wavelengths are slightly above the 
natural frequencies 2ix fn of the cavity while in the attractive case they are slightly below. 



4.2.1 Infinite Barrier 

Let us now consider the limit g — > at fixed n, in which the potential wall becomes infinitely 
high. We find that 

k (n) n, 6k {n) 0, a fc( n) -> 0. (31) 

Then, for the discrete momenta = n, the wavefunction is completely inside the "cavity", 
while for all the remaning values, the wavefunction is completely outside it. The values 
^(n) _ n aDOV e correspond to the quantised (allowed) momenta of a particle in a one- 
dimensional box of length I = n, implying that the potential wall becomes impenetrable in 
the limit g — > 0, in agreement with physical intuition: there is no coupling of the cavity with 
the outside [5]. Outside the cavity, we have instead a continuous spectrum of eigenfunctions 
labelled with k, defined for x > tt and vanishing in x = it. [f| We may say that, in the limit 
g — ^ 0, the system decomposes into two non-interacting sub-systems, representing stable 
particles and a continuum of multi-particle states. It is a remarkable formal fact that in the 
limit g — > the eigenfunctions corresponding to the eigenvalues k^ become normalizable: 
since a fe o) = 0, eq. flTUJ) becomes meaningless and one has to go back to eq. (jSJ) and to 
impose discrete-state normalization — typically normalization to one. 

Let us now consider the properties of the discrete spectrum of our model for g — > 0. For 
g — > + there are no bound states, while for g — > 0~ there is always one bound state which 
gets progressively more localized around the potential support x = tt. Heuristically, this 
discontinuous behaviour of the discrete spectrum around g = suggests that this point a 
singular one. 



5 Time evolution for the unstable state 

This is the central section of the paper, in which we compute the time evolution of unstable 
states by means of non-perturbative analytic techniques. The first problem is therefore that 
of defining an unstable particle and the second one to find the analog of such a state in our 
model. In physical terms, an unstable particle (or resonance) is related to the enhancement 
of any cross section producing that particle in the s-channel in a narrow band of P 2 ~ M 2 , 
where is the total 4-momentum and M the particle mass. An unstable particle is 
therefore a quantum (i.e. coherent) state prepared at t = 0, which is normalizable and 

3 Since k is a continuous variable, removing integer values (a zero measure set) has no influence on the 
spectrum. 
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Figure 1: Integration contour in k-plane for g = 0.1. 



has a non-trivial time evolution, not being an eigenstate of the system @. In general, there 
is some freedom in the definition of such states. We try to keep as close a connection as 
possible with the procedure followed in perturbation theory: we take as unstable state an 
exact eigenstate of the non- interacting system (g = 0), which is no more an exact eigenstate 
of the interacting system (g ^ 0). Let us then consider the evolution of a wave-function 
given at the initial time t = by 

ip m ( x ,t = 0) = J^- sin (I x) 6(tt - x) , (32) 

where I is a positive integer. Note that, unlike the eigenfunctions ipk{%) considered above, 
is a normalized state: 

| tp {l) (x, 0)| 2 dx = 1. (33) 

Let us remark that any continuous wavefunction with support in the interval (0, tt) - 
representing the more general state for an unstable particle — can be expressed as su- 
perposition of the above wavef unctions. As discussed in the previous section, for g = 
the wavefunction above cannot decay being an eigenfunction and therefore we expect it 
to represent a slowly decaying state for < \g\ -C 1. The decay products corresponds to 
non-normalizable wavefunctions in our model of the form 

0fc(x, t) oc 9(x — tt) sin \k(x — tt)\ , (34) 

vanishing for x < tt. For g ^ 0, also the wavefunctions above are not eigenfunctions; in 
physical terms, they can excite modes inside the cavity. The wavefunctions in eq. (|32|) 
exactly vanish outside the barrier, implying there are no "decay products" at t = 0. 
To compute the time evolution, namely 

^ {l) (x, t) = e- i6t ^ l \x, 0), (35) 



4 A general definition of an unstable particle via S-matrix elements constructed by means of wave-packets 
has been provided in [3] (for general background see for example [6]). 
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the most convenient technique is to expand the state above into eigenfunctions of the 
Hamiltonian which, as well known, have trivial time evolution: 

POO 

^ l \x, 0) = / (p®fa(x)dk, (36) 
Jo 

where: 

m® = r-iV; ,, 

satisfying ip_ k = —(f k , with 



<Ph = (-irizNMur- 1 !*' ( 37 ) 



Nk ^ 4a k b k 1 + l/(7igk)sm2k7i + 1/ {2ir 2 g 2 k 2 )(l - cos2£;7r) ' ^ 

We have considered for simplicity's sake only the repulsive case g > 0. Our wavefunction 
at time t is therefore given by: 

POO 

^ l \x, t) = / ^^ k {x)e~ ikH dk. (40) 



5.1 Asymptotic Expansion of Wave Function 

The spectral representation in eigenfunctions of the unstable state at time t has the explicit 
expression: 



2 y/ 2 - 



^ l \x,t)=[-J I p {l \k;x,g)e- tkt dk, < x < vr , g>0, (41) 



where 



sin A/7T ! 

p m (k;x,q) = (-1)H— 77: v, t\ : t, — n n , ow -— sin fcc . (42) 

y V ; K ' k 2 -I 2 l + l/( y ngk)sm2k7r + l/(27i 2 g 2 k 2 )(l-cos2kn) v ' 

To obtain rigorous analytic formulae, we expand the integral for large t. The steepest 
descent method suggests to replace the integral on the r.h.s. of eq. (I4"II) by the integral 
over the steepest descent ray (0, ooe~ l7r ^ 4 ), on which the fast oscillation of the integrand 
is absent. This is achieved by computing the above integral over the sequence of closed 
contours j n containing the segment on the real axis (0, n+a n ), the ray ((n+a n )e~ ln ^ 4: , 0) and 
a circular arc c n connecting the endpoints of the two segments (see fig. 1). The parameter 
a n (0 < a n < 1/2) is chosen is such a way that c n passes at the maximal possible distance 
between the poles in the fourth quadrant of the fc-plane. Therefore the state ip^(x,t) is 
decomposed in a natural way into the sum of two quite different contributions: 

1>®(x,t) = 4 l U^t) + 1>%p(?,t) , (43) 



5 Since the integrand is an even function of k because both ip9 and ip k (x) ar e odd, one can extend the 
integral over all fc's as 

1 f°° 

dk -> - / dk. (39) 
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where 



V&OM) 



9 \ 3/2 



7r y 



-in/4:. 



,x,g) e 



-k 2 t 



dk : 



2 \3/2 oo 

" ' 71=1 



(44) 
(45) 



In general, the contribution ipp6w(x,t) exhibits a power decay as t ^> 1, while the contribu- 
tion ipex P (x,t), coming from the residues in eq. (1431) . exhibits an exponential decay. Let us 
consider the above contributions in turn: 

1. the integral ippow(x,t) is over the ray (0, oo e~ 47r / 4 ) and for large t ^> 1 takes the 
dominant contribution from a neighborhood of k = 0, where the integrand is analytic 
and can therefore be expanded in powers of k: 



P (l) (k;x,g) 



y) i=i 



(46) 



The first few coefficients explicitly read: 

f_iy+i 



pi\x,g) 

pi\x,g) 



I 



Tlx : 



-ixx 



2 ^g 



* 2 g 2 



1 ^ 

P + ~6 ' 31 + g {l + g) 2 



x 



(47) 
(48) 



Replacing this series into the integral and performing the change of variable v = k 2 1, 
one obtains the following asymptotic expansion: 



i> { ±{x,t) 



pj\x,v) 



y/2 e~™! A g 2 



dvv ] 2 e v 



(49) 



V2e^/ 4 g 2 ^ (_i)i(2j - l)!!pf (x,<?) 



vr (1 + ^)2 



E 



2J 



< x < vr , t > 1,(50) 



whose first few terms read: 

e"/ 4 (-1)' g 2 



x 



^2 I (l + g)H^ 



1 



3i 
2t 



1 7T 2 2 7T 2 



/ 2 6 3 1 + # 



X 

~6~ 



(51) 



Let us make a few remarks. The above asymptotic expansion is uniformely valid for 
all g > 0, since the coefficients Pj (x,g) are uniformely bounded in that region (see 
eq. (1461)). The exponent 3/2 controlling the power decay, ip ~ l/t 3//2 , does not depend 
on I and g; 
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2. the simple poles k^ n \g) of the integrand in ipix P (x,t) are the simple zeroes of the 
trascen dental equation 



Aa h b 



sin 2kir 1 — cos 2kir 

1 + — + - ; — = 



(52) 



ngk 2n 2 g 2 k 2 
constrained by the conditions: 

Refc (n) > |ImA; (n) |, Imfc (n) < 0. (53) 

It is easy to prove that all the poles k^ n '(g) satisfy all of the above conditions for 
g <C 1. The second condition in (|5"3"|) is always satisfied. In general, the poles leave 
the fourth quadrant for very large values of \g\, where the unstable-particle description 
is irrelevant. The unstable states (or resonances) have pole contributions for g 1 




Figure 2: Time evolution of the modulus square of the exponential contribution (continuous 
line) and power contribution (dashed line) to the wavefunction of the fundamental state 
I = 1 (see eq.^EM)) for g = 0.2, integrated in the interval (0,7r). The scale on the vertical 
axis is logarithmic. 



of the form: 



where 



with 



1, oo 



CLn = ("I 



d+n 



2ln 



l 2 -n 2 ' 

VftCM) = ^Jlz^(g) sin [k^(g)x]e 



-iu] { - n \g)t-lj2Y {n \g)t 



Z [n \g) 



1-f +0(<7 2 ) 



u {n) (g) = + (Rek^y - (Imk^y 

= n*[l-2g+0(g 2 )} ; 

T {n \g) = -4Refc (n) Imfc (n) 

= 47mV [1-Ag + O (g 2 )} 



(54) 

(55) 
(56) 

(57) 

(58) 
(59) 
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6 Discussion 



Let us consider the asymptotic expansion for t > 1 of the first resonance (I = 1) by 
neglecting the contributions of the higher-order poles (n > 1): 



ip (1 \x,t) 



3 i7r/4 



X 



^2 {l + g) 2 t^ 




7T 2 2 ii 2 g 

1 + — + - 

6 31 



+ 



g 



'2 

7T V 2 

A few comments are in order: 



sin [(1 



x\ 



g (i + g) 2 

exp [-i (1 -2flf)t - 2ng 2 t] . 



x 

~6~ 



1. in the limit g — > 0, the power term disappears from the r.h.s. of eq. (160]) and the 
exponential term approaches the fundamental eigenfunction of a particle with mass 
m = 1/2 in a box of length / = 7r: 



^ (1) (x,t) ->• si 



sin ./• ( . (yf — )• . (61) 

That is in complete agreement with physical intuition, as already discussed; 

2. it is clear that, for sufficiently long times, the decay law will be dominated by the 
power term. However, since the power term has a small coefficient, suppressed as 
g 2 <C 1 for small g, while the exponential term has a coefficient of order one and 
it remains 0(1) as long as t < 1/g 2 , the exponential term dominates over the power 
term for a large temporal region for g ^ 1. In other words, if g is small and t is large, 
but such that 

log(l/<?) 



1 < t < 



(62) 



there is a long transient in which the exponential term prevails on the power term (see 
fig. 2). Since the signal rapidly decays with time, the transient region may actually 
be the one measurable one. 
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NOTE ADDED 

After the first version of this note was put on the archive, references [7] and [8] were brought 
to our attention^. While the first paper deals with general properties of the exponential 
region, the second one treats the time evolution with the saddle point method of the same 
model as we do. We are in complete agreement with [8] as far as the asymptotic power 
behavior in time is concerned, while we are in disagreement with the exponential behavior. 
More specifically, the first term on the r.h.s. of our eq. (I54[) . i.e. the diagonal one n = I, is 

6 We wish to thank Dr. R. Rosenfelder for pointing out to us. 
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in agreement with the r.h.s. of eq.(2a) in [8] (I labels the initial state and n the pole). In [8] 
however, the non-diagonal pole contributions n ^ I, present in eq. (|54p . are not included. 
These terms have a coefficient suppressed by a power of g <C 1 compared to the diagonal 
one, but have a slower exponential decay for n < I (T^ oc n 3 , see eq.(l59l). and therefore 
dominate at intermediate times (i.e. before power-effects take over). As shown in figEJ 
there is indeed a large temporal region where the non-diagonal contribution from the first 
pole, 

Y9 2 £\^ieM\ 2 dx, (63) 
dominates over that of the second pole, 

r \4>SL{x,t)\ 2 dx , (64) 
Jo 

in the temporal evolution of the first excited state, I = 2. Neglecting the non-diagonal 
contributions is therefore a reasonable approximation only for the time-evolution of the 
lowest-lying state 1 = 1. More details will be given in a forthcoming publication [9]. 

0.1 I 



io- 4 - 




50 100 150 200 

Figure 3: Time evolution of the contributions to the I = 2, i.e. first excited, state for g = 
0.1. Dotted line: second pole contribution; Dashed line: first pole contribution; Continuous 
line: power contribution. 
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